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ABSTRACT 

Context. The Herbig Ae/Be star HD 100546 hosts one of the most extensively studied discs in the southern sky. However, not much 

is known about this disc's gas and large scale kinematics. 

Aims. Our aim is to establish whether the disc is gas-rich, and to study the disc temperature and kinematics. 

Methods. We detect and study the molecular gas in the disc at spatial resolution from 7'.'7 to 18.9".' using the Atacama Pathfinder 

Experiment telescope*. The lines '-CO / =7-6, J =6-5, J =3-2, '^^CO J =3-2 and [C I] 'Pa-'Pi are observed, diagnostic of disc 

temperature, size, chemistry, and kinematics. We use simple parametric disc models and a molecular excitation and radiative transfer 

code to fit the observed spectral line profiles. 

Results. Our observations are consistent with more than lO"'' Mq of molecular gas in a disc of r;400 AU radius in Keplerian rotation 

around a 2.5 Mq star, seen at an inclination of 50°. The detected '-CO lines arise from gas at 30-70 K. Asymmetry in the '^CO line 

emission suggests that one side of the outer disc is colder by 10-20 K than the other, and is not due to pointing offsets, cloud absorption 

and asymmetry in the disc extent. 

Conclusions. Both low- and high-y '^CO lines are dominated by the outermost disc regions, and indicate a 400 AU radius. The '-CO 

/ =6-5 line arises from a disc layer higher above disc midplane, and warmer by 15-20 K than the layer emitting the J =3-2 line. The 

existing models of discs around Herbig Ae stars, asuming a B9.5 type model stellar atmosphere overproduce the [CI] 'P2-'Pi line 

intensity from HD 100546 by an order of magnitude. 

Key words. Planetary systems:protoplanetary disks - stars: individual: (HD 100546) - stars: pre-main-sequence 



W^ 1 . Introduction ture dWilner et alJ I2003L iRaman et al.l l2006t iPanic et al.l l2008l 

O l2009h . 
^^ Over the past decade our understanding of the structural and 

^ physical properties of discs around young stars has increased ^ntd recently, observations of rotational transitions of 

^: from basic theoretical modelling of the spectral energy dis- molecules in the submillimetre regime were focused primar- 

•^ tributions (SEDs) constrained by observations with no spatial ^^^ "" *« l"^"-^ emission from '^CQ, up to the 7 ^3-2 line 

X information, to modelUng based on not only the SEDs, but (Greaves et al. 2000; Thi et al. | | 200U i Qi et alj | 2004 |Thi et al. | 

H also spatially resolved dust observations, like scattered light ^QQ4^i Dent et al. ,| 2005aD . In two of the brightest and most stud- 

. P. . images and interfe rometry (Pinte et a l. 2008; P anic et al.llSt '^^ ^o™'^^^' ^W Hya and LkCa 15, the observations of higher-/ 

iTannirkulametai] l2008l) . Two decades ago, the first submil- transitions of '^CO, up to J=6-5 (E^ =116 K), were compared 

Hmetre interferometer observations resolved the molecular gas ^o the low-/ lines, providing estimates of the gas temperature 

emission spatially and this allowed major progress in under- mjhe mtermediate-height molecular layer (| van Zadelhoff et al. | 

standing the disc k inemati cs, structure a nd che mistry (e.g., l2Q01a|), crucial ingredients for chemical modelling of discs. 

Beckwith & Sarg^i^ flQSTt iKoerner et al.1 flool iDutrev et al.1 ^^^^ single-dish line spectra were fitted using simphstic disc 



1994D.Amorere^t example is a study of one of t he brightest '^°^^^''^ ^^^"^^"8 ^ temperature of 20-40 K in the '^CO lineemit 



discs around a Herbig Ae star, HD 163296 shown in llsellaet al.1 t/IlKiayers of LkCa 15, and more than 40 K in TW Hya.i^ielalJ 

(l2007h . (Sub)millimetre gas and dust emission is the ideal probe ^^W^ analysed submillimetnc interferometer observations of 

of the global disc properties, like size, mass and radial distri- ^W Hya m the context of a disc structure based on an accre- 

bution of disc material, because the bulk of the disc mass is ^i"" ^isc model ^alvet et al. | |2002| ). Based on '^CO J =6-5, 

located beyond 100 AU from the star, at temperatures of 10- ^-2 and 2-1 observations, they show that X-ray heating of the 

50 K that dominate this part of the specti-um. Disc models 8^^ is efficient in this source, in addition to the stellar radiation 

which include consti-aints of both dust and molecular gas ob- ^^1^. Such diagnostics of gas heating and lonisation improve our 

servations have sti-essed the importance of analysing the gas and understanding of how the gas content evolves in discs, 
dust components simultaneously, in the context of a single disc The emerging (sub)millimetre facilities in the Southern 

model with a three-dimensional temperature and density struc- hemisphere like the Atacama Pathfinder Experiment (APEX) 
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and the Australia Telescope Compact Array (ATCA) are open- 
ing a window towards the star-forming regions of the Southern 
sky and are well suited to study circumstellar disc emission. 
These instruments also pave the path for future observations 
with the Atacama Large Millimetre / Submillimeter Array 
(ALMA), which will drastically improv e our knowledge of 
disc structure and evolution (e.g. iGuillot eau & Dutrev 2008). 
We use APEX receivers APEX-2a and CHAMP+ to observe 
the '^CO J =7-6, J =6-5, J =3-2, '^co J =3-2 and [C 
I] ^P2-^Pi line emission towards the disc around the young 
intermediate-mass star HD 100546. A wealth of observations 
of dtist in this bright di sc (WaeUcens et al. 1996; Malfait et al. 
19981: iGradv et al.ll200U IXu gereau et al. 2001; Bouwman et al. 



2003 



Acke & van den Anckeri i2006ai: lArdila et al.i i200 7) has 



motivated us to probe its molecular gas content and kinemat- 
ics. The chosen transitions are particularly sensitive to the gas 
in the warm upper layers and kinematics of the outer disc. 
Our millimetre line observations probe the outer radius and in- 
clination. The existing observational constraints on these pa- 
rameters in the disc around HD 100546 provide an excellent 
basis for the analysis of our data. Our observations also pro- 
vide a bridge toward even higher-/ far infrared '^CO tran- 
sitions to be observed with the Herschel Space Observatory. 
HD 100546 is a young B9V type, 2.5 Mq stai-, classified as 
a Herbig Be star due to its isolation, infrared excess and sili- 
cate emission (iThe et al.lll994l; iMalfait et al.l ll998). With a dis- 
tance of 103+6 pc, measured by Hipparcos, this is one of the 
nearest Herbig Ae/Be stars. The age of the star is est imated 
to be greater than 10 Myr (va n den Ancker et al.l Il998 i). This 
makes the presence of circumstellar material intriguing, consid- 
ering that discs are found to dissipate within 10 Myr in most 
young stars (e.g..l Hollenbach et al.ll2000t [Hernan dez et a l.|2007l: 
[Hillenbrand 2008). Based on SED modelling, Bouwm an et al.l 
(1200 3) postulate the presence of an inner hole in the disk with 
a 10 AU radius, possibly c aused by a Jupiter-sized planet (see 
also lAcke & van den Ancke r 2006a). Direct evidence of cold 
disc material at larger radii is provided by ATCA observations of 
IWilneret all (12003) at 89 GHz (3.4 mm) and 2" resolution, with 
the flux of 36+3 mJy, v alues consistent with the 1.3 mm observa- 
tions of lHenning e"tan(il998 ). They do not detect HCO+ 7=1-0 
line emission and speculate that photodissociation of its progen- 
itor species ^^CO, in the upper disc layers or an overall gas de- 
pletion may be the reason for this. Recent spect roastrometric ob- 
servations of ro vibrational '^CO transitions by v an der Plas et alj 
(.2009.) suggest that this molecular species is missin g from the in- 
ner di sc, at least up to 8 AU from the star (also see lBrittain et alJ 
l2009h . 

Scattered light imaging of HD 100546 reveals the disc ex- 
tending up to 4" from the star viewed at an inclination of 
50°, and an intere s ting disc structure resembling spiral arm s 
dPantin et all 120001; lAugereau et all 120011: iGradvetalJ l200lh . 
This structure yvas interpreted as d ue to disc perturbation by 
a co mpanion ( Quiflen et al.l l2005h or a war ped disc struc- 
ture (Ouillenll2006t) . Coronographic imaging bv lAugereauet alj 
(1200 ih shows steep surface brightness profiles in the environ- 
ment of HD 100546 indicative of optically thin emission in the 
near-infrared, with surface densities as low as 10"^ g cm"^. Their 
images trace the emission of small dust (< 5jt/m), extending out 
to 800 AU from the star The authors suggest the presence of an 
optically thick disc with a 400 AU radius, and an optically thin 
flattened halo or envelope farther from the star The scattered 
light images hint at the presence of gas in the disc that supports 
the disc vertical structure. 



In this work, we detect and study the molecular gas, its kine- 
matics and temperature, in the disc at spatial resolution from 
7'.'7 to 18.9'.'. In Sect. |2] we present our observations of '^CO, 
'^CO and [C I] lines. All '^CO lines are detected, there is a ten- 
tative detection of the '^CO line, while [C I] emission is not 
detected. We model the spectra in Sect. |3] discussing the impli- 
cations for the disc size, mass and kinematics. We identify the 
regions dominating the observed lines, and derive their temper- 
atures. Section|4]summarises our results. 

2. Observations and results 

The observations of '^CO J =6-5 at 69 1 .472 GHz and [C I] ^P2- 
''Pi ([C I] 7 =2-1 hereafter) at 809.344 GHz towards HD 100546 
at R.A.=11*'33"25!4 and Dec= -70°11'41" (J2000) were ob- 
tained simultaneously, u sing the CHAMP^ heterodyne array re- 
ceiver dGusten et al.l200 6) on APEX on 2008 November 1 1 . The 
7 pixels in each wavelength band are arranged in a hexagon of 6 
pixels around one central pixel pointed towards the source, with 
beam sizes of 9" at 691 GHz and 7'.'7 at 809 GHz. The data were 
obtained in the position switching mode in which the telescope 
moved between the on-source position and a reference position 
at an equatorial offset of (-1000", 0") from the source. The 
backend consisting of the Fast Fourier Transform Spectrometer 
unit was used on all pixels, providing a spectral resolution of 
0.37 MHz or 0.14-0.16 km s"' at these frequencies, and cover- 
ing a bandwidth of 1.5 GHz or 4096 channels. Main beam effi- 
ciencies are 0.40 at 691 GHz and 0.37 at 809 GHz. Calibration is 
uncertain by s;30% at both frequencies. Pointing was performed 
directly prior to, and after the on-source integration providing an 
accuracy better than 3". The pointing source was 07454-7112, 
about 20° away from HD 100546, at the same airmass. The CO 
J =6-5 line was also observed on 2008 November 10 in jig- 
gle mode and its intensity and spectral profile were found to be 
the same within 20%. Since these data are noisier, the data sets 
taken on different dates were not combined. During this latter 
observation, the high band was tuned to '^CO J =7-6 at 806.665 
GHz, with the beam efficiency of 0.37. Only the high S/N '^CO 
J =6-5 data taken on 2008 November 1 1 are used in the further 
analysis. 

The I^CO J =3-2 line at 345.796 GHz, and the '^CO J =3- 
2 line at 330.588 GHz, were observed on 2005 July 27 and 
28 with the APEX-2a receiver using a single pointing. The 
channel spacing of these data is 61 kHz or 0.05-0. 06 km s~', 
and th e spectral resolution 98 kHz or 0.09 km s ' ( jKlein et al.l 
I2006h . The main beam efficiency of APEX at 346 GHz is 0.73, 
and the beam sizes are 18'.'1 and 18'.'9 at 345.796 GHz and 
330.588 GHz, respectively. Calibration is uncertain by a;20% at 
th ese frequencies. Our ' ^CO J =3-2 line data were presented 
in ' Panic & Hogerheiidd (l2009l) . In all our observations, APEX 
forward efficiency of 0.95 was taken into account. 

The data were reduced and analysed using the CLASS and 
STARLINK software. We detect all the observed '^CO line tran- 
sitions at 7-20cr in terms of ffie integrated intensity. The '^CO 
J=3-2 line is marginally detected at 2cr, and the [C I] J =2-1 line 
is not detected. Figure [T] shows the observed spectra, corrected 
for the beam efficiency and re-binned to a lower spectral resolu- 
tion. A first degree polynomial fit, except for the 0-10 km/s range 
in which the lines are emitted, is subtracted from the full spec- 
tral range for each data set. The '^CO J =3-2 and J =6-5 lines 
are detected at the highest signal to noise ratio, and show a dou- 
ble peaked profile characteristic of disc rotation. The intensities 
integrated over the velocity range 0-10 km s"', over which line 
emission is detected, are listed in Table [1] together with the full 
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Table 1. Integrated intensities of the observed spectral lines. 



Spectral Line 


/" 


FWHM 




(Kkms-') 


(km s-i) 


''CO J =1-6 


12.9±1.9 


4.2 


''CO J =6-5 


17.7±0.9 


4.2 


"CO J =3-2 


4.0±0.6 


4.0 


'^CO J =3-2 


1.3+0.6* 


(...) 


CI J =2-1 


<1.0^ 


(...) 



Notes. 

*"' The line intensity / = / Tub dV is integrated over the velocity range 

0-10 km s'. 

**' Marginal 2cr detection. 

*''' Upper limit given by the 2cr value. 



width at half-maximum (FWHM) of the lines with sufficiently 
well defined profiles. 

In addition to the observations towards the source, the 
CHAMP^ array provides measurements at nearby offsets. This 
setup provides an excellent way to discern the emission from the 
disc, with an estimated size of at least 400 AU in radius from 
the surr ounding material known to be present further away from 
the st ar dHenninget alJll998t iGradv et al]l200lL lAugereau eTaP 
1200 ih . The 400 AU lower limit to th e size of the gas disc is 
set by the scattered light observations JAugereau et al.ll2001h is 
close to the spatial resolution of our data. The central pixel of our 
CHAMP^ data probes the '^CO J =6-5 line emission from the 
region of 9" centered on the position of the star (460 AU radius), 
while the surrounding pixels probe the more distant regions (at 
a distance 1000-2000 AU from the star). Similarly, at the fre- 
quency of the '^CO J =7-6 line a smaller region around the star 
of 7'.'7 is probed with the central pixel (400 AU radius), and re- 
gions roughly 1000-2000 AU from the star with the surrounding 
pixels. Table |2] provides an overview of the pixel positions and 
the corresponding fluxes integrated over a velocity range from 
to 10 km s ' velocity range, over which the '^CO lines are 
firmly detected in the central pixel. Compared to the on-source 
fluxes, these measurements show that the '^CO emission from 
the surrounding material is about 20 times weaker than from the 
region within 400 AU from the star. The '^CO J =3-2 line was 
observed with a single pointing, and the beam of IS'.'l is large 
enough to include any emission from regions beyond 400 AU. 
However, the strong resemblance in the line profile suggests that 
both low-7 and high-7 lines arise from the disc and that any con- 
tribution to the line emission by an extended low-temperature 
and low-velocity component is negligible. 

Our data testify to the presence of warm molecular gas at 
distances less than roughly 400-500 AU from the star, corre- 
sponding to the smallest beam sizes in our observations (7'.'7 
and 9"). We obtain the following integrated intensity ratios, cor- 
rected for beam dilution (scaled to the same beam): '^CO lines 
(6-5)/(3-2)=l.l±0.6 and ('^CO 3-2)/(i^^CO 3-2)>2.8+2.2. The 
eiTors include both the rms and flux calibration uncertainty. 

There is a clear asymmetry in the profile of the '^CO J =3-2 
and J =6-5 lines, observed at a high signal-to-noise ratio, char- 
acterised by a 20-40% difference in the line intensity between 
the two peaks. This asymmetry is certain, given that the spectra 
shown in Fig. [T]have been smoothed to a lower spectral reso- 
lution, and that each spectral bin contains an average of several 
measurements. 

In some sources, asymmetries like this are explained through 
confusion with foregro und cloud emission, as seen in '^CO 
lines towards IM Lup dvan Kempen et all l2007h and DL Tau 




-20-15-10-5 5 10 15 20 25 30 

V 



Vlsr ('<m s ) 



Fig. 1. From top to bottom: Spectra of the '^CO J =7-6, J =6- 
5, J =3-2, '^CO J =3-2 and CI J =2-1 lines observed towards 
HD 100546 shown with black lines plotted as a histogram. The 
'^CO and '^CO spectra computed from our models of the line 
emitting layers of the disc (see Sections 13.2.11 and 13.2.21 ) are 
shown with continuous red lines. At 100 AU from the star, our 
vertically isothermal models of the '^CO line emission have tem- 
peratures rioo=60-70 K, while for the '^CO a lower temperature 
7^100=30 K is used. 



dSimon et alJl2000h . However, in our off-pixels of the CHAMP+ 
array we see no significant extended emission. While a slight 
pointing offset in the 6-5 line observations towards a symmet- 
ric source may result in a line asymmetry, the persistance of 
line asymmetry in our IS'.'l resolution observations of the 3-2 
line indicates that the asymmetry is due to the source properties. 
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Table 2. The '^CO J =6-5 and J -1-6 integrated line intensities at offset positions. 



'-CO J =6-5 


'^CO J =7-6 


RA Offset 


Dec Offset 


^CO(6-5) 


RA Offset 


Dec Offset 


-^00(7-6) 


(") 


(") 


(Kkms-') 


(") 


(") 


(Kkms-') 


-9.8 


-17.0 


0.4 ± 0.2 


+ 10.2 


-17.0 


<2.0 


+8.8 


-18.0 


0.5 ± 0.2 


+20.2 


-1.4 


<2.0 


-18.9 


-0.2 


<0.2 


-9.0 


-16.5 


<2.0 


-H19.4 


-1.2 


<0.2 


+ 11.0 


+ 16.2 


<2.0 


-9.3 


+ 16.6 


<0.2 


-18.8 


+0.2 


<2.0 


+ 10.0 


+ 16.6 


<0.2 


-9.1 


+ 16.9 


<2.0 



Notes. The line intensities / = J ^mb dV are integrated over a velocity range 0-10 km s ' . The upper limits are given by the 2cr value. 



Furthermore, the two observations were taken years apart and 
a systematic offset in both observations is highly unlikely. We 
therefore exclude pointing offset as a cause of the observed line 
asymmetry. If one side of the disc extends slightly further out 
than the other (e.g., on the SE side), the increase in disc surface at 
that side will contribute to the line flux at one side of the spectral 
line causing line asymmetry. This may be a plausible explana- 
tion for the 3-2 line, but not the observed 6-5 line with a 9" beam 
size, unless the asymmetry is within 450 AU. Such an asymme- 
try may affect the scattered light images of HD 100546, with 
the disc extending further to the SE than to the NW of the star 
Howev er this is not seen in the observations of lAugereau et alJ 
(1200 ll) . Furthermore, the additional emission from this region 
would dominate the frequencies closer to the line centre, and not 
the blue-shifted peak. 

Asymmetry is seen in '^CO lines from sources with a pro- 
nounce d structural disc asymmetry, like HD 141569 dDent et alj 
l2005bl) . Another possibility is a temperature asymmetry in the 
disc. In the following section we analyse and model the observed 
lines and investigate these scenarios. 

3. Discussion 

3.1. Line ratios as temperature and density diagnostics 

3.1 .1 . '^CO 3-2/'3CO 3-2 line ratio and disc density 

The submillimetre '^CO line emission from circumstellar discs 
is generally optically thick and arises from warm upper disc lay- 
ers. These lines are particularly sensitive to the temperature of 
these layers, and therefore to the stellar and external illumina- 
tion of the disc surface. Comparison to isotopologue line emis- 
sion is often used to assess the optical thickness. If the ^^CO 
3-2 emission is treated as a detection, we obtain the line ra- 
tio of (i^CO 3-2)/('3CO 3-2)=2.8+2.2. We make a rough es- 
timate of the column density needed to reproduce the derive d 
ratio using the RADEX online toolQ (I van der Tak et all 120071) . 
For kinetic temperatures 20-40 K, H2 densities 10^-10^ c m~^, 
and an isotopic ratio {^^CW^C\^11 (Wi lson & Roodll 1994 . the 
('^CO 3-2)/('3CO 3-2) ratio derived above indicates the ^^CO 
column density of roughly 10'^ cm"^, and excludes lower col- 
umn densities. Using a standard '^CO abundance of 10 ■* with 
respect to H2, this translates to the gas column density larger 
than 10^^ cmr^. 

3.1 .2. '^CO e-S/'^CO 3-2 line ratio and disc temperature 

We compare the line ratio (^^CO b-5)lC^CO 3-2)= 1.1 +0.6 to 
the calculations done using RADEX, and assume that these two 



' |http : //www ■ strw . leidenuni v . nl/~Bioldata/radex . html | 



lines are dominated by the same region in the disc. This assump- 
tion is generally t rue up to 200-300 AU fro m the star (see cal- 
culations done in Ivan Zadelhoff et al.ll2001b.) . For H2 densities 
of 10''-10'^ cm -3, and '^CO column densities of lO'^-lO'"* cm ^ 
the ('^co 6-5)/('^CO 3-2) line ratio a;l.l is reproduced at tem- 
peratures of 60-100 K. If the disc is large and sufficiently cold, 
the assumption that the two lines arise from the same region is 
not necessarily true, as the 3-2 line may be dominated by colder 
material at radii larger than the 6-5 line. It is therefore uncer- 
tain whether the line ratio reflects the difference in the temper- 
ature between the two vertical layers or in the temperature of 
the radial regions emitting the lines. We explore this further in 
Sect. 13. 2.TI using more sophisticated radiative transfer tools and 
parametric models of disc structure. 

The observed CO 6-5/3-2 integrated intensity ratio of 
1.1+0.6 (observations scaled by the beam size) is higher than 
the ratios close to 0.5 found for t he discs around the TTauri 
stars LkCa 15 and TW Hya by Ivan ZadelhoffeFai] ( 1200 lal) 
and lOi et alJ (l2006l) . The TW Hya 6-5/3-2 ratio has been 
interpreted as proof that the gas temperature is higher than that 
of the dust in the surface layers where gas and dust are not 
thermally coupled. Both UV radiation and X-r ays have been 
invoked to provide the additional gas heating ("Jonkheid et alj 
^004; Kamp & Duflemond 2004; Glassgold et al. 200i 
iNomura & MilJail l2005t iGorti & HoUenbachI 12008). although 
for H erbig Ae stars the X-rays can be neglected (.Kamp et al] 
I2OO8I) . 

The higher ratio found for HD 100546 implies higher gas 
temperatures than for the aforementioned disks around T Tauri 
stars. This is expected based on models where most of the gas 
heating comes from UV radiation from th e central star, sinc e 
the cooler stars have less UV radiation (e.g.. lWoitke et alJ2009l) . 
The gas temperature is also strongly affected by the PAH abun- 
dance in the disc. R4H features are seen pr ominently in th e 
HD 100546 mid-infrared spectrum (e.g., Ma lfait et al.lll998l) . 
and have been shown to be spatially extended across the disk 
dOeers et al. 2007). In contrast, neither of the two T Tauri disks 
for which the '-CO J =6-5 line has been detected show PAH 
emissi on, implying typic al PAH abundances a factor of 10-100 
lower dGeersetalfcooel) . 

Fo r the specific case of Herbig Ae disks, iJonkheid et al.l 
(I2OO7I) have computed the gas heating and chemistry as well as 
the resulting CO 6-5 and 3-2 li ne emission starting from a set o f 
dust disk models developed bv iDuIIemond & DominikI (l2005h . 
Models with decreasing disk mass from 10"' to 10"'' M© and 
decreasing dust/gas ratios from 10"^ to 10"^ (simulating grain 
growth and settling) were investigated. The UV radiation field 
was t aken to be that of a B9.5 star (see Fig. 3 of Ijonkheid et al.l 
I2OO6I) . very close to that expected for HD 100546. An accu- 
rate treatment of the shape of the UV field at short wavelengths. 
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<1 100 A, is very important for a correct calculation of the CO 
photodissociation and atomic carbon photoionization rates, since 
the UV intensity of a B9 star is orders of magnitude weaker in 
this wavelength range than the (scaled) standard interstellar ra- 
diation field. The PAH abundance is taken to follow the dust/gas 
ratio, with an abundance of 10"^ for a normal dust/gas ratio of 
0.01. 

The ^^CO J -6-5/3 -2 integrated intensity ratios computed 
bv lJonkheid et alJ (l2007h summed over the full extent of the disk 
model are remarkably close to unity for the entire range of disk 
parameters investigated, consistent with our observations. The 
absolute '^CO and '^CO intensities for model B2 (a disk with 
a mass of 0.01 Mq with a standard gas/dust ratio of 100, ap- 
propriate for HD 100546) are also within 40% of the observed 
values when scaled to the same source distance and beam size, 
indicating a good agreement between models and observations. 

3.2. Parametric disc models and origin of ^^CO line emission 

Submillimetre '^CO line emission is analysed using two differ- 
ent modelling approaches in the literature. Dis c physical mod- 
els, lik e the irradi ated accretion di sc mo dels of iD'Alessio et alj 
(l2005h : Dullemo nd & Dominid (l2005h . are especially well 
suited when the emission is spatially resolved. In this way, 
the three-dimensional structure of the disc can be investi- 
gated, for exa mple when transit ions of different optical depths 
are observed dPanic et al.l l2008h . For spatially unresolved ob- 
servations, like those presented here, simplistic models with 
a limited number of free parameters are more appropriate 
to derive some basic constraints on disc propert ies based on 
the l i ne spectrum (method emp l oyed in, e.g.. 
19941 



. , ^ iDutrev et al] 

Guillo teau & Dutrev* '1998*: 'v an Zadelhoff et alj 1200131 



Dartois et al.l 2003; Dent et al. 2005 a). Here we go a step fur- 
ther than in Sect. l3.1l in analysing the observations in the context 
of such models. 

De nsity and temperature structure. iPanic & Hogerheiidd 
(l2009h show that simple power-law disc models, with a disc 
mass M - 0.01 M©, an H2 surface density S oc /? ' and tem- 
perature T = Tioo K(/?/100 AU)"« with Tioo =60 K and q =0.5, 
are a useful tool to analyse low-7 '^CO transitions from discs 
around gas-rich Herbig Ae stars. These models are vertically 
isothermal and have a vertical density structure in hydrostatic 
equilibrium. The temperature and surface density of such model 
with an outer radius of 400 AU is shown in Fig.|2] upper panels. 
The '^CO abundance with respect to H2 is assumed to be 10""*, 
constant throughout the disc. Their analysis uses the assumption 
that the low-/ '^CO lines arise from the upper disc layers (opti- 
cally thick) and that their fluxes are dominated by the outermost 
disc regions. Therefore the spectral profiles of these lines are 
very sensitive to the extent of the emission, with strong and nar- 
row lines indicating a large disc and weak and relatively wider 
double-peaked lines characteristic of small discs. 

Spatial origin of the emission. The regions of the disc 
probed by the ^'^CO / =6~ 5 and J =3-2 lines are studied in 
Ivan Zadelhoff et al] (1200 Ibl) . Their Fig. 6 shows the locations of 
the surfaces of opacity t=1 for these lines in a cross-section of 
their disc model. At distances up to 150 AU from the star, the 
two transitions trace similar depths into the discs, but at larger 
radii they trace different vertical layers: the low-excitation 3-2 
line probes deeper and colder layers than the 6-5 line. To as- 
sess whether both 6-5 and 3-2 line fluxes are dominated by the 
same radial disc region, and to identify this region, we use the 
calculations of synthetic 6-5 and 3-2 spectra from disc mod- 
els as described in the previous paragraph. We use two differ- 



ent temperature profiles, with q =0.5 and q =0.2. The radial 
temperature and surface density structure of the two models is 
shown in Fig. |2l In this exercise we use the face-on disc orien- 
tation, for simplicity, and calculate the emerging '^CO J =6- 
5 and J -li-l line intensity along lines of sight at a number 
of different radial locations from the star This is done using 
the molecular excitation an d radiative transfer code RATRAN 
(iHogerheiide & van der Tak 2000) and CO-H2 collision rates 
from the Leiden At omic and Molecular Databas^ (LAMBDA; 
ISchoier et al.l 120051) . We are primarily interested in knowing if 
there is a single radial region dominating the line flux more than 
other regions or the flux is distributed more evenly with the ra- 
dius. Also, we investigate how the regions of strongest line flux 
spatially compare for the two '^CO transitions. 

The contribution to the line flux from an annulus of a fixed 
width dR located at a distance R from the star is given by 
dF{R) = 2nI{R)RdR, and therefore dF{R) cc I{R)R. The rela- 
tive contribution of the annuli to the total line flux arising from 
the model dF{R)/F,o,ai is also proportional to I{R) R, allowing us 
to directly compare the importance of one radial region with re- 
spect to another. The Figure|2]shows dF{R) for the '^CO J =6-5 
(dashed line) and J =3-2 (dotted line) in case q =0.5 (upper 
panels) and q =0.2 (lower panels). For ease of comparison be- 
tween the two lines, the '"CO J =3-2 and J =6-5 line fluxes 
are plotted with respect to their maximum value dF^^^. For both 
temperature profiles, the outermost disc region (200-400 AU) 
is very important, especially for the 3-2 line. For this line the 
strongest flux arises at the largest radii due to the efficient ex- 
citation at low temperatures found in the outer disc combined 
with the larg e surface area. This result is in line with our previ- 
ous results in lPanic & Hogerheiidd ( l2009h where the '^CO 7=3- 
2 spectra are used as probes of disc sizes. While the 6-5 line 
is also sensitive to the outer disc regions, the upper panels of 
Fig. |2] show that our model assumptions at 50 AU matter just 
as much as those at 400 AU. This is explained by the excita- 
tion energy of the '^CO 7=6-5 line, that is somewhat high for 
the 30-40 K temperatures found in the outer disc in the model 
with a steeper temperature profile ^=0.5. The flatter tempera- 
ture distribution with T >40 K in the outer disc allows to ef- 
ficiently populate the higher energy levels of '^CO, and in this 
case the outermost disc regions are clearly more important than 
50-100 AU, due to the mere increase of the surface area with ra- 
dius. Understanding the diagnostic potential of the high-7 '^CO 
lines in discs is useful already for the analysis of spatially unre- 
solved APEX observations (see Sect. 13.2.1b . but also shows that 
the future sub-arcsecond observations of high-7 '^CO lines in 
discs with ALMA could be used to probe the radial temperature 
structure in discs. 

Disc mass. Our initial choice of parameters, described above, 
roughly matches the observed '^CO J =6-5 and J =3-2 line 
fluxes, confirming the outer radius of 400 AU seen in scattered 
light observations. We performed test calculations for a lower 
value of the disc mass of lO""' Mq, i.e., ten times lower den- 
sities. The modelled line fluxes are found not to vary signifi- 
cantly with respect to the noise levels of our observations. Thus, 
in the disc mass range of lO^^'-lO"^ Mq, values found typically 
towards circumstellar discs, and temperatures similar to those 
assumed here, the observed '^CO J =6-5 and J =3-2 lines are 
insensitive to disc mass. If we decrease the disc mass to allow 
the '^CO line emission in our models to be sensitive to the den- 
sity variations in the midplane (i.e., emission dominated by the 
midplane), the disc midplane temperatures (<20 K in the outer 
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Fig. 2. (a): Radial temperature structure of a disc model from lPanic & Hogerheiidd (2009), with /?out =400 AU, Tioo =60 K and 
q =0.5 is shown with the black full line. The radial dependence of the line flux dF(R)/dF(R)„^^^ for the '^CO 3-2 line and 6-5 
line arising from this model is overplotted using the dashed and the dotted lines, respectively.ffe).- The surface density of the same 
model is shown with the black full line, and line fluxes as in (a), (c) and (d): Same as in (a) and (b), for the model with a shallower 
temperature profile, q =0.2 and its corresponding line fluxes as a function of radius. 



disc) would be the only physically appropriate assumption for 
the temperature in the 'low-mass' model. While the asymmetry 
in the line profile could be reproduced in this way, the resulting 
lines would be much weaker than observed. Given these consid- 
erations, we assume the two lines to be optically thick in the disc 
around HD 100546 and fix the disc mass to the initial value of 
10"^ Mq. Our model is consistent with the thermal continuum 
emission of d ust observed t oward HD 100546 by Henning et al. 
(IT998); Wiln er et all ( l2003h . In the further analysis we fit the ex- 
act line shapes and intensities, with a particular focus on the line 
asymmetry seen in both J =6-5 and J =3-2 spectra. 

3.2.1 . i^CO J =6-5 and J =3-2 line fit 

The characteristic temperatures and masses derived in Sect. 13.11 
fo r the disc a round HD 10054 6 are similar to the models used 
in lPanic & Ho gerheiid d (l2009h. The modell i ng of the '^CO 3-2 
line of HD 100546 in Panic & Hogerheiidd (l2009l) yielded very 
rough estimates of the outer radius of 300 AU and inclination 
of 35°, consistent with the 400 AU and 50° from the scattered 
light. We use their models to fit the '^CO spectra, with Tioo as 
free parameter and varying q from its initial value of 0.5 where 



necessary. We fix the outer radius and inclinatio n to the observa- 
tionally constrained values 400 AU and 50° (see lAugereau et al] 
2001, and the discussion of disc size in Sect. |2]l. We use the 
Keplerian velocity field around a 2.5 M© star The dependence 
of the spectral profile on the stellar mass in the range 2.0-3.0 Mq 
does not affect our fit significantly. 

The inner radius is assumed to be 0.6 AU, close to the dust 
sublimation radius. Although an inner hole of 13 AU is found 
in HD 100546, its presence would not aff'ect our results because 
the molecular lines observed are dominated by disc regions far 
beyond the inner tens of AU. Because the observed '^CO lines 
trace warm molecular material (>20 K) and are insensitive to the 
colder regions deeper in the disc, it is reasonable to neglect the 
'^CO freeze-out. 

As above, the line emission is calculated using RATRAN. 
The synthetic images from our disc models are convolved with 
the beam sizes of the corresponding observations, and the spec- 
tra are extracted towards the centre of the image. Where tem- 
perature asymmetry is modelled (see 13.2.1b . pairs of synthetic 
images with diff'erent Tioo parameter were first combined - con- 
tributing each to a given side of the disc with respect to the semi- 
major axis. Following this, the combined image is convolved and 
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spectra extracted. Dust continuum emission, although negligible 
for the molecular line transfer, is included in the calculation and 
subtracted from the synthetic image cubes. 

In the view of the result that the line emission appears to be 
insensitive to the density, the observed line asymmetry can not 
be caused by a density asymmetry in the disc, but it rather re- 
flects an asymmetry in the disc temperature structure. To model 
such an asymmetry, we use the fact that the emission at the two 
peaks in the line profile is dominated by the two sides of the disc 
with respect to the minor axis. We assume that one side of the 
disc is colder than the other. This may be due to a combined ef- 
fect of an asymmetry in disc illumination and geometry. We use 
different Tioo parameters for the two sides of the disc and we 
explore different radial temperature distributions with T oc /? ?, 
^=0-0.5. 

The best fit to the '^CO J =3-2 spectrum, seen in Fig. [1] 
is obtained by assuming the radia l temp erature structure with 
^=0.5, as in Panic & Hoge rheiidd (|2009|) . and Tioo of 60 and 
70 K for the two disc sides, respectively. The line shape is con- 
sistent with the assumed Keplerian velocity field and disc incli- 
nation. The temperatures of 60-70 K compare well to the theoret- 
ical predictions of temperatur e in the regions where thes e lines 
saturate in discs (see Fig. 6 in Ivan Zadelhoff et al.ll2001ah . They 
show that at 100 AU from the star the J =6-5 and / =3-2 lines 
trace the same disc layer, we assume that the '^CO J =6-5 line 
is emitted at Tioo of 60 and 70 K, as found for the 3-2 line. We 
find the best fit to the '^CO J =6-5 spectrum, shown in Fig. [1] 
assuming T =60 Kx(/?/100AC/)""-' on one side of the disc and 
T =70 Kx(/?/100 At/)-" 2 Qjj jhe other Thus, the radial temper- 
ature distribution in the 6-5 line emitting la yer is flatter than in 
the 3-2 layer, consistent with the result of Ivan Zadelhoff et alJ 
(1200 Ibl) for the outer disc regions where the 6-5 line traces 
warmer material higher in the disc surface than the 3-2 line. 

Therefore our '^CO line observations, although spatially un- 
resolved, provide some information on the temperature radial 
distribution in the 6-5 and 3-2 line emitting layers and disc tem- 
perature asymmetry. It is important to note that slightly different 
sets of temperature values and slopes than those found here may 
provide an equally good fit to these observations. However, any 
such model requires higher outer disc temperatures to fit the 6-5 
line than required for the 3-2 line, clearly showing that the two 
lines arise from different vertical layers in the outermost disc re- 
gions. 

The difl^erence in temperature between the two sides of the 
disc may be explained by a warped inner disc, as illustrated in 
Fig. [3] with the elevated side of the inner disc intercepting a frac- 
tion of stellar light that would otherwise reach the outer disc, 
while the opposite side is slightly more illuminated, with the in- 
ner part of the disc shift ed downwards. The possibility of an in- 
ner warp is suggested in lOuilleril (l2006l) for HD 100546, with an 
inner component extending up to 200 AU inclined by a; 15° with 
respect to the outer component extending beyond that radius. 

The temperature asymmetry is also possible if the disc has a 
different thickness at the two sides. This may happen in a disc 
with dust settling underway, if a planet or another body embed- 
ded in the disc stirrs the dust back up. In this case the 'stirred' 
part of the disc intercepts more stellar light and becomes some- 
what warmer A companion body in HD 100546 is suggested in 
the literature, to e xplain the observed i n ner hole, gas kinematics, 
and sp iral arms ( lOuillen et al.l 120051: lAcke & van den Anckeij 
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Fig. 3. Illustration of a warped disc structure, in an edge-on pro- 
jection. Straight light and dark grey lines show the plane of the 
disc and of the inner warped disc component, respectively. The 
white curve shows the surface region with a higher temperature, 
heated well by the star, and the black curve shows the cooler re- 
gion partly obscured by the warp. The image of the disc with a 
400 AU radius is adapted from Ouillea (.2006.) . where a detailed 
warped disc model is presented. 



3.2.2. i^co J =3-2 and '^CO J =7-6 line fit 

With the low signal to noise in our '^CO J =3-2 spectrum 
(marginal detection at 2a-), it is not possible to assess the spec- 
tral line shape. The disc temperatures used to fit the '^CO lines 
in Sect. 13.2.1! overestimate the '^CO J =3-2 line intensity. A 
lower temperature of the '^CO J =3-2 emitting layers with re- 
spect to the '^CO J =3-2 (as discussed in Sect. 13.2.2b . a higher 
['^C]/['-^C] isotopic ratio than the standard interstellar medium 
value of 77, or an overall decrease of CO abundance can explain 
our marginal detection of the '^CO line. Freeze-out and/or se- 
lective photodissociation processes in discs are known to affect 
the CO abundance and isotopic ratios. Detailed modelling and 
spatially resolved submillimetre line observations of '^CO and 
isotopologues would allow to constrain the disc structure, eval- 
uate the effect of these processes (,e.g.. P anic et al. 2008) and 
constrain the amount of molecular gas better. Clearly, there are 
a number of uncertainties involved in the interpretation of our 
'^CO J =3-2 data. For these reasons we limit the analysis to a 
comparison of our marginal detection to one simple model. A 
lower temperature of the '■'CO emitting layers is expected due to 
a different optical depth of '^CO lines with respect to the '^CO. 
In Fig. [T] we compare our data to model calculations done as- 
suming low temperatures, rioo=30 K and q-Q.5, an isotopic ra- 
tio ['^C]/['-'C]=77, and all remaining model parameters as in the 
simple axially symmetric power-law model shown in Fig. |2] and 
described in Sect. l3.2.T] Considering that the '^^CO line emission 
probes colder and deeper disc layers, it is likely that temperature 
asymmetry derived for the upper layers from '^CO lines is neg- 
ligible for '^CO. High-sensitivity observations of '^CO or C'^0 
submillimetre lines in the future may reveal symmetric line pro- 
files. 
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The '^CO J =7-6 and J =6-5 transitions arise from energy 
levels at /s^ of 155 K and 1 16 K, respectively, and are expected 
to trace roughly the same disc layers. Figure [1] shows how our 
best-fit model of the 6-5 line compares to the observed 7-6 line. 
The observed line flux is about 50% lower than the model pre- 
diction. This may be due to a slight mispointing of the telescope, 
as the size of the source (800 AU in diameter) is close to the 7'.'7 
beam size at this frequency. The noise levels of these observa- 
tions are notably higher than in the 6-5 line observations and 
spectral profile not well defined. Due to these reasons, we do 
not fit the 7-6 line further but use the more reliable 6-5 line fit 
results in our analysis. 

3.3. Implications of the [C I] J -2-1 non-detection 

Figure 1 includes the high-quality spectrum around the [C I] 
J =2-1 line at 809.344 GHz. No significant feature is detected 
down to 0.3 K rms in a 0.27 km s ' velocity bin, implying a 
2 cr limit on the integrated intensity of » 1.0 K km s"' over 
the range 0-10 km s~' (same width as for the detected '^CO 
lines). Model B2 of lJonkheid et al.l(l2007h predicts integrated [C 
I] intensities, scaled to the distance of HD 100546, around 15- 
20 K km s ', whereas the model intensities are even larger in 
disks with significant grain growth and settling (BL model se- 
ries in the abovementioned paper). Thus, while the CO data ap- 
pear entirely consistent with their sophisticated UV-heated disk 
atmosphere models, the [C I] data are clearly discrepant by an 
order of magnitude. 

One possible solution could be that the radiation field con- 
tains more carbon ionising photons than assumed here, shift- 
ing the chemical balance from neutral to ionised atomic car- 
bon. The predicted [C II] line intensities for the model disks of 
[jonkheid et al. ( 2007) are very low, <0. 1 K km s"' , whereas they 
are more than an order of magnitude hi gher for disks aroun d 
T Tauri stars with excess UV emission (.Jonkheid et al.l 120041) . 
consistent with an eflicient carbon ionisation around Herbig Ae 
stars. Such excess of UV emission could come from the disk- 
star accretion boundary layer. Indeed, HD 100546 is observed to 
undergo significa nt accretion, in spite of the known (dust) gap 
in the inner disk (IVieira et al.l[T99 9). Also, the FUSE UV spec- 
trum of HD 100546 shows significant flux, both in the contin- 
uum and in lines such as O VI, in the wavelength range where 
atomic carbon can be ionised (Lecavelier des Etangs et al. 2003, 
and G. Herczeg, private communication). The continuum emis- 
sion should also photodissociate CO, but not the stellar UV lines 
unless they overlap with the discrete CO photodissociating tran- 
sitions. Also, because CO is self-shielding and the lines are op- 
tically thick, the CO emission may be less affected than that 
of atomic carbon by any extra UV. However, detailed chemical 
modelling is required to test this hypothesis. 

[C II] 158 i-im emission has been detected by Malfait et al. 
(1998) with ISO-LWS and most recendy with the PACS instru- 
ment on the Herschel Space Observatory by Sturm et al. (2010, 
submitted). The ISO fluxes are up to an order of magnitude 
larger than the model fluxes by Jonkheid et al. (2007) but a sig- 
nificant fraction of the emission observed in the large ISO beam 
may be caused by general Galactic background emission. More 
detailed modeling is needed to pin down the origin of the ab- 
sence of the [C I] 809 fim line and check whether the observed 
[C II] flux is consistent with enhanced UV radiation. Future ob- 
servations of [C I] J =l-01ineat492.161 GHz (^Pi-3Po), which 
is generally brighter than the [C I] J =2-1 line, may provide a 
stronger constraint on the amount of neutral atomic carbon in 
the disc. Observations of the [C II] line with Herschel-HIFI can 



directly resolve the line profile and determine whether this line 
indeed comes from the disc. Overall, the [C I] flux and [C I]/[C 
II] line ratio could become interesting diagnostics of the UV ra- 
diation field to which the disc is exposed. 



4. Conclusions 

We summarise our conclusions as follows: 

- We present evidence for warm molecular gas associated with 
the disc around HD 100546, in the regions within 400 AU 
from the star, successfully separated from more extended 
material in our CHAMP^ observations; 

- The gas kinematics are consistent with Keplerian rotation 
around an 2.5 Mq star of a disc with a 400 AU radius, viewed 
at an inclination of 50° from face-on; 

- The I^CO (6-5)/(3-2) line ratio of 1.1+0.6 is higher than 
measured towards discs around T Tauri stars, likely due to 
a more efficient gas heating of the disc containing PAHs by 
the stronger UV radiation from the B9 star; 

- Our data testify to the significant molecular gas reservoir in 
the disc, consistent with the total disc masses of more than 
lO""* Mq. We exclude the possibility of a low-density disc 
and optically thin '^CO emission. 

- Line asymmetry seen in '^CO J =6-5 and J =3-2 lines is 
explained by a temperature asymmetry, with one side of the 
disc slightly colder than the other, possibly due to a partial 
obscuration of one side by a warped inner disc or a high disc 
rim. We exclude radial asymmetry, midplane density asym- 
metry and mispointing as possible causes; 

- Our modelling shows that, due to the efficient heating of 
the disc gas by the star, both low- J and high-7 '^CO lines 
are dominated by the outermost regions of the disc, though 
slightly different vertical disc layers with Ar= 15-20 K; 

- We find that in 'colder' discs where temperatures of the emit- 
ting regions are close to 20-30 K in the outer disc - colder 
Herbig Ae discs and especially T Tauri discs - the high-/ 
lines probe a larger extent of the disc, starting from as little 
as 50 AU; 

- The non detection of [C I] J =2-1 line may indicate the pres- 
ence of more carbon-ionising photons than assumed in the 
B9 model atmosphere. 

Future observations with ALMA will be crucial to charac- 
terise the disc around HD 100546, and spatially resolve its kine- 
matics and structure. In particular, these observations will al- 
low a detailed comparison between the spatial distribution of the 
gas traced by the rotational transitions of '^CO and its isotopo- 
logues, and the dust traced with the millimetre continuum emis- 
sion. Herschel far-infrared data can probe even higher J '^CO 
transitions, as well as [O I] and [C II] lines originating from the 
disc surface. Being a bright, isolated source suspected to harbour 
a planet, the disc around HD 100546 is one of the prime targets 
to probe disc structure in the early planet-building phases. 
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